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Abstract 

We calculate the cross sections for neutrino interactions with (isoscalar) nuclear targets 
in the energy domain all the way up to 10 12 GeV. Small x QCD effects are included by us- 
ing a unified BFKL/DGLAP formalism which embodies non-leading log 1/x contributions. 
The few free parameters which specify the input parton distributions are determined by 
fitting to HERA deep inelastic data. The attenuation of neutrinos transversing the Earth 
at different nadir angles is calculated for a variety of energy spectra for neutrinos origi- 
nating from different sources (from Active Galactic Nuclei, Gamma ray bursts, top-down 
models), as well as for atmospheric neutrinos. For this purpose we solve the transport 
equation which includes regeneration due to neutral current neutrino interactions, besides 
attenuation. 



1. Introduction 

The penetration of ultrahigh energy neutrinos through the Earth, with energies E greater 
than 10 TeV or so, can be strongly affected by neutrino interactions with matter. This is due 
to the increase of neutrino cross sections with energy. At these ultrahigh energies we have 
significant attenuation of the neutrino fluxes on transversing through the Earth and, indeed, 
complete absorption for energies above about 10 8 GeV or so, depending on the nadir angle of 
the neutrino beam. Realistic estimates of these effects are crucial for predicting the number of 
neutrinos reaching the large km 3 scale detectors after penetration through the Earth. Clearly 
the magnitude of this ultrahigh energy neutrino flux has important implications for neutrino 
astronomy. For example, it is hoped that the neutrino flux coming through the Earth will point 
back to its Active Galactic Nuclei (AGN) origin |IJ. Clearly the magnitude of the signal, and 
the ability of the Earth to reduce the background due to atmospheric neutrinos, are crucial in 
this exciting endeavour. 

The inelastic interaction of neutrinos with nucleons is traditionally described by the QCD- 
improved parton model. In ultrahigh energy neutrino interactions we are probing a kinematical 
region which is not accessible in current collider experiments. To be precise the most powerful 
electron-proton collider HERA at the DESY laboratory in Hamburg typically probes the region 
x > 10~ 4 for Q 2 > 10 GeV 2 , whereas the ultrahigh energy neutrino interactions can become 
sensitive to the domain x ~ 10~ 8 and Q 2 ~ M^. Here, as usual, x denotes the Bjorken 
scaling variable and Q 2 = —q 2 , where q is the four momentum transfer between the leptons 
in the inelastic scattering of a lepton on a nucleon. The Bjorken variable x is defined as 
x = Q 2 1 {2p.q) where p denotes the four momentum of the nucleon. 

In perturbative QCD it is expected that the gluon and sea-quark distributions, and hence 
also the deep inelastic scattering structure functions, should grow with the decreasing values of 
the parameter x. This theoretical expectation has been beautifully confirmed by the structure 
function measurements at HERA . These measurements put important constraints on parton 
distributions in the small x region probed at HERA. In order to get predictions for the ultrahigh 
energy neutrino-nucleon cross sections one has to construct a reliable extrapolation of the 
structure functions to the region of very small values of x which is probed in these ultrahigh 
energy interactions. The existing extrapolations, which do also incorporate the constraints 
from the HERA data, are based entirely on leading order (LO) or next-to-leading order (NLO) 
DGLAP evolution [^|, [|, [|. This approximation may, however, be incomplete at low x since it 
ignores the important resummation of ln(l/x) BFKL-type effects |§. In this paper we wish to 
incorporate the QCD expectations at low x which will take these effects into account. To be 
precise we shall base our calculation of the neutrino cross sections on the unified BFKL/DGLAP 
formalism which incorporates both the ln(l/x) resummation and the complete LO DGLAP 
evolution J7|. Due to the size of the NLO ln(l/x) contributions || one might question the 
reliability of this procedure. However our framework makes it possible to resum dominant 
non-leading ln(l/x) contributions to all orders. This has the important effect of stabilizing the 
non-leading contribution and turns out to give a physically and phenomenologically acceptable 
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suppression of the LO BFKL behaviour. In this way we should achieve the most reliable 
dynamically-motivated extrapolation to very low x to date, which incorporates all the relevant 
QCD expectations. The calculated neutrino cross sections are then used as an input in transport 
equation describing the propagation of the neutrinos through the Earth R TCj. This equation 
will contain both the attenuation of the neutrino flux together with its regeneration through 
the neutral current interaction. Similiar calculation which takes into account absorption and 



gives prediction for the muon rates has been performed in [11 



The content of our paper is as follows. In the next Section we discuss deep inelastic lepton- 
nucleon scattering at low x within the unified BFKL/DGLAP framework. In Sec. 3 we collect 
the relevant formulas needed for calculating the neutrino cross sections and present our nu- 
merical results for these cross sections calculated within the unified BFKL/DGLAP scheme. 
We also confront our predictions with results of calculations based on NLO DGLAP evolution. 
Sec. 4 is devoted to the discussion of the transport equation and its solution for the shadowing 
factor due to passage of the neutrinos through the Earth. In Sec. 5 we present our results 
concerning the changes of the initial fluxes with nadir angle. We consider neutrino fluxes cor- 
responding to Active Galactic Nuclei, Gamma Ray Bursts and the Top-Down models, together 
with atmospheric neutrino background. Finally in Sec. 6 we give our conclusions. 



2. Deep inelastic scattering at low x 

Ultrahigh energy neutrino-nucleon interactions probe values of Bjorken x which can be 
several orders of magnitude smaller than those which are accessible at the present deep inelastic 
ep scattering experiments at HERA. Here, as usual, x = Q 2 /2p.q where Q 2 = —q 2 , with p 
denoting the nucleon 4-momentum and q being the 4-momentum transfer between the leptons 
in the deep- inelastic process £N — > £'X. 

At low values of x we must consider log(l/x) effects. In ref. |7j a formalism is presented 
which permits an extrapolation of parton distributions to very small x. Besides the leading 
order (LO) log(l/x) resummation, the procedure incorporates a major part of an all-order 
resummation. We outline the method below. We begin with the BFKL equation || for the 
unintegrated gluon distribution f(x,k^), which performs the LO log(l/x) resummation. The 
equation has the form 

/ex,*) = /<»>(*, *»> + * s e r±f d ^S f^p-{ir /z ' k2) + ID 



x Z 



k' 2 [ \k' 2 -k 2 \ [4k' 4 + k A ]^ 



where as — N c as/^ and k = kr,k' = k' T denote the transverse momenta of the gluons, see 
Fig. 1. The term in the integrand containing f(x/z,k' 2 ) corresponds to real gluon emission, 
whereas the terms involving f(x/z,k 2 ) represent the virtual contributions and lead to the 
Reggeization of the t-channel exchanged gluons. The inhomogeneous driving term will be 
specified later. 
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The observable nucleon structure functions Fi are given in terms of the gluon distribution 
by the kx factorization formula [fjj| 

F^Q 2 ) = £ ^ / ^i^ box (.,fc 2 ,g 2 )/(^ 2 ), (2) 

where F^ ox describes the subprocess Vg — > qq, see Fig. 1. Here the virtual gauge boson V 
may be either a 7 (describing an electromagnetic deep inelastic scattering) or a boson 
(describing a charged current weak interaction) or a Z boson (describing a neutral current 
weak interaction). The procedure (P and (||) automatically resums all the leading log(l/x) 
contributions to the observable Fi. 

The solution of the LO BFKL equation for fixed as gives a QCD or hard pomeron with 
intercept a(0) = 1 + A with A = cF,s41n2. The ln(l/x) resummation has recently been carried 
out m at next-to- leading order (NLO). It is found to give a very large 0(a s ) correction to A 

A ~ a 5 4 In 2(1 -6a s ), (3) 

which implies that the NLO approximation is unreliable for realistic values of as- Rather we 
must use a formalism which contains an estimate of an all-order resummation. Clearly it would 
be desirable to identify physical effects which could be resummed to all orders and which at 
the same time yield a NLO value of A that is comparable to (13). As it happens the imposition 



of the consistency constraint [12, 13 



k' 2 < k 2 /z (4) 

on the real gluon emission term gives just such an effect. The variables are shown in Fig. 1. 
The origin of the constraint is the requirement that the virtuality of the exchanged gluon is 
dominated by its transverse momentum \k' 2 \ ~ k^. For clarity we have restored the subscript 
T in this equation. 

If condition (^) is imposed on the BFKL equation it can be still solved analytically. The 
result is an all-order effect, which at NLO gives the large modification 

A ~ a 5 4 In 2(1 - 4.2a s ) (5) 

of the LO value. However it is found that the all-order correction is a much milder modification, 
although still significant. A related result can be found in ref. [H]. We can therefore make the 
BFKL equation (0) for the gluon much more realistic by imposing the consistency condition 
(f|), as well as by allowing the coupling as to run. 

Moreover we can extend its validity to cover the full range of x. To do this we note that the 
BFKL equation embodies the important double leading log part of DGLAP evolution which is 
driven just by the singular 1/z part of the splitting function P gg . To obtain a reliable description 
throughout the full x range (and not just at small x) we must include the remaining terms in 
Pgg, together with the quark to gluon transitions. We also introduce in eq. [I] the parameter k 2 
(k 2 ~ lGeV 2 ) which divides the non-perturbative ikl 2 < k 2 ) from the perturbative (k' 2 > k 2 ) 
region. Finally we note that the contribution from the infrared region k' 2 < k% in ([!]) may be 
expressed [0] in terms of the integrated gluon distribution at scale k 2 ,, that is g(x, fcg). All the 
above modifications of (0) are encapsulated in a unified BFKL/DGLAP equation of the form 
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f(x,k 2 ) = f {0 \x,k 2 ) + 



J x Z J k 



dk 



rl 



f[-,k' 2 

z 



,12 



kl k' 2 



\k' 2 - k 2 \ 



+ — — ~T 



[Ak' A + k A }\ 



(6) 



rX dz ( z 



PM - 1 



+ as(k 2 ) [ 
Now the driving term has the form 



2tt 



dzP gq (z)T, 



(0), 



2vr 



dzP gg (z)-g \-,kl 

z \z 



(7) 



where in the perturbative domain k 2 > k 2 , we may safely neglect the genuinely non-perturbative 
contribution f^°'(x, k 2 ) which is expected to decrease rapidly with increasing k 2 . It is important 
to note that (^|) only involves f(x,k 2 ) in the perturbative region k 2 > k%. The input (|7|) is 
provided by the conventional gluon at scale &q, just as in pure DGLAP evolution. 

The last term in (Q) is the contribution of the singlet quark distribution to the gluon, with 



E = $>(9 + 9) = E(^ + K 



(8) 



where S and V denote the sea and valence quark momentum distributions. The gluon, in turn, 
helps to drive the sea quark distribution through the g — > qq transition. Thus equation (§) has 
to be solved simultaneously with an equivalent equation for E(x, k 2 ). 

Consider for the moment just the g —>■ qq contribution to S q . The hp factorization theorem 
gives [|15[ 

S q (x,Q 2 ) = ^ *Z j ^Sl x (z,k 2 ,Q 2 )f(- z ,k 2 ) (9) 

where S hox describes the quark box (and crossed-box) contribution shown in Fig. 1. S 
implicitly includes an integration over the transverse momentum k of the exchanged quark. We 
have 



S* ox (z, k 2 , Q 2 



Q 2 



4ir 2 k 2 



d(3 / dVa s < [(3 2 + (1 - f3) 2 ] 



K 



k — k 



D 



D 



2q 



+ [ml + 4Q 2 (3 2 (l-(3) 2 } 



D 



D 



2<Z, 



6(z- z ) 



(10) 
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where k' = K — (1 — f3)k and 

D lq 



= k 2 + /3(1-/?)Q 2 + m\ 



D 2q = {n-kf + (3{l-(3)Q 2 + m 2 



20 



/t' 2 + m\ k 2 
+ (3{l-(3)Q 2 + Q~ 2 



(11) 



Eqs. (P|)-(P^|) enable us to evaluate the singlet quark distribution £ in terms of the gluon /. 
We obtain 

£(*, k 2 ) = S„(x) + E T - S ^ k ' 2 = °> fc2 )- 9 (-, k 2 

„ J x Z Z \ Z 

r°° dk' 2 r l dz 



(12) 



dk' 2 a s (k' 2 ) f^^ fx u2 



where a = (1 + Am q /Q~) and V = x(u v +d v ). Here we have separated off the non-perturbative 



contributions. S^n-p is the contribution from the region k , k < k 2 and the next term is 
the contribution from the region k 2 < k 2 , < k,' 2 . The details are explained in ref. [[/J. An 
S — > S contribution (from the light u, d, s quarks) is also included. For the light u, d, s quarks 
Sq OX (z, k' 2 = 0, k 2 ) in (|12|) is defined with the n' 2 integration restricted to the region k! 2 > k 2 . 

In this way coupled integral equations are obtained for the unintegrated gluon f(x, k 2 ) 
and the integrated quark singlet S(x, k 2 ) distributions. The driving terms are specified by 
an economically parametrized non-perturbative contribution S non ^ p (x) and by the integrated 
gluon distribution at scale = 1 GeV 2 which was taken to be of the form 

xg{x,kl) = N(l-xf. (13) 

The valence distribution V(x, k 2 ) was taken from the parton set of ref. |l6fl . With this input 



the coupled equations (@) and ([12]) were solved and a fit made of the deep inelastic electron- 
proton F2 data obtained by the HI, ZEUS, NMC and BCDMS collaborations 0. An excellent 
description of these data was obtained with physically reasonable values of the parameters: 
= 1.57 and (3 = 2.5. Incidentally, the output gluon at x ~ 0.4 was reasonably compatible 
with the expectations of prompt photon data. In summary, ref. gives as unintegrated gluon 
distribution and, through kx factorization, a sea distribution which can be reliably extrapolated 
to very low values of x. 
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3. The neutrino cross sections 

Here we collect together all the relevant formulas which are needed to calculate the cross sec- 
tions for neutrino-nucleon interactions. We express the cross sections in terms of the structure 
functions for an isoscalar nucleon target, N = (n + p)/2, |17|, [18], [19| 



d 2 ^ G F ME ( M[ \ f l + (l-y) 2 2 , 
dxdy vr \Q 2 + M 2 \ 2 j 



V l F l{x,Q 2 )±v [!-"-) xF»(x,Q 2 



(14) 



where Gp is the Fermi coupling constant, M is the proton mass, E is the laboratory energy of 
the neutrino and y = Q 2 /xs. The mass Mj is either or Mz according to whether we are 
calculating charged current (CC) or neutral current (NC) neutrino interactions. 

We may express the structure functions Fi in terms of the valence and sea quark momentum 
distributions, V q and S q , of (|T0|) . For the charged current vN interaction we have 

F 2 CC = x(u v + d v ) + S u + S d + 2x(s + c + b + t), (15) 

where the heavy quark contributions are calculated from the photon-gluon fusion mechanism. 
To be precise for the s —> c and c — > s contributions we use 

2*g(a;, Q 2 ) = r (fc2=0) - S^(z, k 2 = 0, Q 2 ) * g k 2 ) 

J x Z Z \ Z J 

(16) 



g k 2 Jx z q v ' ' J \z 



with q = s or c and with Sg OX defined by ( |T0D with m q = 0. However the mass is included in 
the threshold factor 



a c (k 



2\ 



k 2 + m 21 1 

1 + 



(17) 



Q 2 

For the small b — > t and i — > b contributions we use the standard on-shell factorization formula 

2xq(x,Q 2 ) = j — H(z } m q ,m q t,Q 2 )- g(-,s) (18) 

Jx Z Z \Z J 

with q = b or t (and q' — t or 6) and where the scale s = Q 2 (l — z)/z and 

a = [l + ( mt + mb )7Q 2 ]-i. (19) 

The functions H are defined in ref. 
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(20) 



The function B q describes the boson-gluon fusion contribution to Fl 

of a s(Q 2 ) f 1 dy x ( x\ 2 

B q {x,Q) = / [1- -} yg(y,k ) 

Jx y y V V 



(21) 



where D iq and zq are given in ([□]). Besides B u and B d , there is also boson-gluon production 
of heavy quarks which may be treated similarly to that for -F 2 CC . 

The structure functions for the neutral current vN interaction are 



if c = (Ll + L 2 d + R 2 u + R 2 d )^J2S q + xu v + xd v }, 
where the sum over S q is given by the first three terms on the right hand side of QT2]), and 



(22) 



if c = (Ll + L d -Rl-R 2 d )\{u v + d v 



(23) 

(Ll + L 2 d + Rl + R 2 d )-J2 B q + [ dy Q " F^(y, Q 2 ) 

where the chiral couplings are 



7 NC 



L u = 1 - | sin W , L d = -1 + | sin 9 W , 



R u = — | sin 2 6y/, Rd = I sin 2 ^^- 



(24) 



The parton distributions that are used to evaluate the neutrino cross sections are those 
described in Section 2, which are obtained by solving the unified BFKL/DGLAP equations with 
higher order In 1/x effects incorporated via the consistency condition. The various components 
of the CC and NC neutrino cross sections are shown as a function of the laboratory neutrino 
energy in Figs. 2 and 3 respectively. In each case we see that for neutrino energies E > 10 5 GeV 
the sea quark contributions dominate over the valence. The rise of the sea distributions with 
decreasing x is reflected in the continued rise of the sea components of the cross section with 
energy. We also note that for E > 10 5 GeV the valence component of the cross section becomes 
independent of energy. This results from the fact that when (^) is integrated over x and y, or 



to be precise over x and Q 2 , the Q 2 integration is effectively cut-off at Q 2 ~ M 2 , together with 
the fact that the number of valence quarks is finite. The threshold effects of the heavy quark 
contributions are also evident in Figs. 2 and 3. We combine the charged current and neutral 
current cross sections in Fig. 4. 

In ref. we also solved the unified equation with the omission of non-leading effects arising 
from the consistency condition. Using those partons we obtain higher cross sections at ultrahigh 



7 



energies as illustrated by the dashed curve in Fig. 5. Of course this is only shown for comparison 
since it is based on LO BFKL which is an unreliable approximation M. Therefore throughout 
this paper we impose the consistency constraint (f|) which generates the dominant non-leading 
ln(l/a;) effects. In Fig. 6 a comparison is made with other recent calculations of the vN total 
cross section @, P based on NLO DGLAP evolution, with BFKL and higher order ln(l/x) 
effects neglected. We see that our results and those of [f|, are remarkably similar considering 
their different dynamic origin. Since we include resummation of \n(l/x) effects we expect 
the continuous curve to give the most reliable extrapolation to ultrahigh energies. In fact 
we may conclude from a comparison of (^) and ([5]) that our all-order sub-leading hi(l/x) 
treatment gives an upper estimate of the cross sections. Considering that the predictions are so 
similar, although they are based on different dynamical assumptions, we may conclude that the 
ambiguity in extrapolating the neutrino cross sections to ultrahigh energies is less than might 
at first be expected. 

At ultrahigh energies the antineutrino cross sections are essentially identical to the neutrino 
cross sections, since the difference is due to the structure function F3 which is controlled by 
valence quarks. In Table I we list the various cross sections as a function of energy. 

To gain insight into the (x, Q 2 ) domain that is sampled by ultrahigh energy neutrinos we 
show plots in Fig. 7 in which the contours enclose regions contributing different fractions of 
the total cross section. The two plots illustrate the dependence of the domains on the neutrino 
energy. For the ultrahigh energy chosen for Fig. 7(b) we see that the main contribution comes 
from the domain Q 2 ~ and x ~ M^ V /(2ME), where M is the proton mass, as anticipated 



from fll4l). For the lower energy used in Fig. 7(a) some residual propagator effects are still 



apparent and Q 2 < and x < M$ V /2ME. 

4. Transport equations 

The ultrahigh energy neutrinos when penetrating through the Earth can undergo attenu- 
ation due to charged and neutral current interactions as well as the regeneration due to the 
neutral current interactions at higher energies. Both effects are summarized in the transport 
equation for the neutrino flux I(E, r) || |TU| : 

= -a m (E)I(E,r) + I JO- ^f ' »> /(g , T ) (25) 
get J 1 — y ay 

where otot = ^cc + °nc an d where y is, as usual, the fractional energy loss such that 

Ef=^- v . (26) 

The variable r is the number density of nucleons n integrated along a path of length z through 
the Earth 

t = f* dz'n(z'). (27) 



o 
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The number density n(z) is defined as n(z) = N&p(z) where p(z) is the density of Earth along 
the neutrino path length z and Na is the Avogadro number. Clearly the number of nucleons 
r encountered along the path z depends upon the nadir angle 9 between the normal to the 
Earth's surface (passing through the detector) and the direction of the neutrino beam incident 
on the detector. For example 6 = 0° corresponds to a beam transversing the diameter of the 
Earth. To compute the variation of r with the angle 9 we need to know the density profile of 



the Earth. We use the preliminary Earth model [21 



It is convenient to represent the solution of the transport equation ( |25| ) in the form: 

I(E, r) = I Q {E) exp(-a TOT (E)r) V(E, r) (28) 

where Io(E) denotes the flux of neutrinos incident on the surface of the Earth from outer space. 
The function fy(E,r) would be unity in the absence of regeneration. It satisfies the following 
equation: 

dV{E,r) f dy I (E') da^ c {E',y) 

j = ~, — — - exp[—(a TOT (E ) — <j tot (E))t\ V(E ,r), 29 

dr J 1 y I (E) dy 

with the initial condition ty(z,r = 0) = 1. We solve this equation numerically and determine 
the regeneration factor \1/ as a function of E and r. The solution is sensitive to the behaviour 
of otot for E' > E and on energy dependence of the initial flux Iq{E) and on the value of 
dcr^c /dy. In particular the flatter the initial spectrum Iq(E) the more it is possible to sample 
da/dy at energies E' much higher than E and so the amount of regeneration is increased. In 
practice the result is a combined effect of the fall-off due to Iq and the experimental attenuation 
and the enhancement due to da^c/dy. To illustrate the general properties of the solution we 
show the shadowing factor, 

S = *{E,t) ex P (-a TOT (E)r), (30) 

of (|28| ) for two different, but physically relevant, forms of the initial flux. First we consider the 
flux of atmospheric neutrinos given by 



I (E) = cE~ 3 - 6 , (31) 
which has a rapid fall-off with energy, and second, we consider the flux from Active Galactic 



Nuclei (AGN) as given by ref. for which the incident flux Io(E) is approximately constant 
throughout the interval 10 3 < E < 10 5 GeV. The results are presented by the continuous 
curves in Fig. 8 for three different incident directions, 9 = 0°,40° and 80°, corresponding to 
values of t/Na = 1.1 x 10 10 , 0.45 x 10 10 , and 0.072 x 10 10 cmwe respectively. For illustration we 
also show, by dashed curves, the pure attenuation factors A = exp(— o"tot t ) with regeneration 
omitted. Since the neutrino cross sections increase with energy we observe that the attenuation 
factor A leads to total shadowing once the energy is sufficiently high. The energy at which 
this occurs depends mainly on r. The regeneration, which increases the flux at the detector, is 
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sensitive to the energy dependence of I . For a steeply falling flux, corresponding, for example, 
to atmospheric neutrinos, regeneration gives a rather small effect. On the other hand for a 
flatter initial flux Iq(E) the regeneration effect can be quite significant. In fact it can even 
enhance the initial flux by as much as 40%, as can be seen from the shadowing factor at 9 = 0° 
for the AGN flux used in Fig. 6. A similar result was found in ref. f9|. In the next section we 
present the flux arriving at the detector from various sources for a range of angles taking into 
account the full shadowing factor of (p8|) . 



5. Penetration of the Earth for given incident fluxes 

The initial incident neutrino flux Iq(E) is modified on its passage through the Earth by the 
shadowing factor S of (|30|) . For experimental purposes the relevant quantity is the neutrino flux 
1(E) reaching the detector at different nadir angles 9. The four plots in each of Figs. 9 and 10 
show first the initial flux Iq(E) and then the flux at the detector 1(E) for the three nadir angles 
9 = 80°, 40° and 0°. Recall that 0° is penetration of neutrinos through the centre of the Earth. 
The various curves are for neutrinos of different origin. For convenience of reference we show 
the atmospheric neutrino flux in both Figs. 9 and 10. In Fig. 9 we compare this background 
spectrum with three different models of the flux expected from Active Galactic Nuclei (AGN) 
23j |24] , [25H . The AGN flux stands out above the atmospheric neutrino background for neutrino 



energies above about 10 5 GeV. However the AGN spectrum is attenuated at ultrahigh neutrino 
energies by shadowing. The smaller the nadir angle 9 the greater the shadowing. For example 
at 9 = 80 or 40° the weighted flux EdN/dE falls below 10" 15 cm" 2 s" 1 sr" 1 at E ~ 10 s and 
10 7 GeV respectively. Fig. 10 shows the corresponding fluxes for neutrinos coming from gamma 
ray bursts |26| and from a sample top-down model P7j. Similar attenuation can be observed 
in Fig. 10 to that in Fig. 9. 

So far we have discussed the effects generated by inelastic neutrino interactions with hadronic 
matter. Another possible source of modification of the neutrino fluxes which penetrate the Earth 
may be due to neutrino oscillations p8|, [29|, BO, [HI] . The properties of neutrino oscillations in 



matter have been discussed in [3^j. They are characterised by an effective oscillation length 
£ m and mixing angle 9 m which differ from those in vacuum, which we denote £ v and 9 V . To be 
precise in the MSW model for mixing of two neutrino species we have 



tan 29 ri 



£ £ 2 
1 + 2 f cos9 v + ^ 

sin 29 v 



cos29 v +(£ v /£) 

where £ originates from the matter contribution to the oscillation length and is given by 

v/2tt 1.77 x 10 7 



(32) 



(33) 



G F n e 



m 



Pt 



(34) 
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where n e is the electron density and p e is the number density in units of Avogadro's number/cm 3 . 
An important property of £ m is that it saturates at the value £ given by ([34]) p3| . That is 

(Umax = I ^ 10 3 - 10 4 km (35) 

for matter oscillations in the Earth with p e = 2 — 10. This is in contrast with the oscillation 
length in vacuum 

£ v = AnE/Aml (36) 

which increases with increasing energy. Most importantly we see that at sufficiently high energy 
£ v /£ becomes very large and the effective mixing angle 9 m tends to zero, see d33|). In fact for 
the very small values of the mass difference of the neutrinos (Am^ ~ 1CT 6 eV 2 suggested by the 
data [0) we can safely neglect the modification of the neutrino fluxes due to neutrino oscilla- 
tions in their passage through the Earth for E > 1 TeV or so. Of course neutrino oscillations 
can, in principle, modify the flux Io(E) arriving at the Earth. 

6. Conclusions 

In this paper we have extended the recently developed []7| unified BFKL/DGLAP framework 
to calculate the total cross sections of ultrahigh neutrino interactions with nucleons. The 
framework incorporated non-leading ln(l/a;) effects which are generated by the consistency 
condition given in (^). In this way it was possible to resum to all orders the dominant part of 
the non-leading ln(l/x) effects and to obtain a physically acceptable description of the structure 
functions at low x. Indeed this unified BFKL/DGLAP approach gave an excellent description 
of the F2 data from HERA and enables us to extrapolate the structure functions to the region 
of the very small values of x which are probed in the ultrahigh energy neutrino interactions. We 
compared our predictions for the cross sections with the results of the two recent calculations 
which were obtained within the NLO DGLAP framework. We find that all three approaches 
give results for the neutrino cross sections within 30% - 40% or so. We may conclude that the 
potential ambiguities in the extrapolation of the cross sections are much smaller than might 
have been expected. However the present calculation, which includes a treatment of ln(l/x) 
effects at small x, should be the more reliable for ultrahigh energy neutrino interactions. It 
should be stressed that the inclusion of the dominant non-leading ln(l/x) effects was crucial 
for obtaining this result. Extrapolation based on the LO BFKL equation would generate cross 
sections which would be enhanced by a factor of more than two at ultrahigh neutrino energies 
E v ~ 10 l 2GeV . This LO approximation is, however, known to be unreliable. 

The calculated neutrino cross sections were next used as an input in the transport equation 
describing the modification of the neutrino flux during the penetration of the Earth. This 
equation incorporated both the attenuation effects of the neutrino "beam" as well as the regen- 
eration of neutrinos due to neutral current interactions. We solved the equation for a variety 
of initial neutrino fluxes and discussed the dependence of the shadowing factor upon the nadir 
angle 9. We found that although attenuation of the neutrino flux is the major effect, neverthe- 
less enhancement due to neutrino regeneration can become appreciable for fluxes which extend 
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to the very high neutrino energies, like those originating from AGN sources. The calculations 
of the neutrino fluxes at various nadir angles show that at sufficiently high energies neutrinos 
become strongly attenuated. The smaller the nadir angle the lower the energy of complete 
attenuation. 

To sum up we have demonstrated that a framework which incorporates QCD expectations 
at low x, including the BFKL effects with the resummation of the non-leading ln(l/a;) terms, 
gives neutrino cross sections which are compatible with those obtained within the NLO DGLAP 
framework. This strongly limits potential ambiguities in the possible values of the cross sec- 
tions extrapolated from the HERA domain to the region of x and Q 2 which can be probed in 
ultrahigh energy neutrino interactions. Due to large values of these cross sections the attenu- 
ation effects reduce the fluxes of ultrahigh energy neutrinos particularly at small nadir angles. 
Nevertheless there is a window for the observation of AGN by km 3 underground detectors of 
the energetic decay muons. We have found that the AGN flux exceeds the atmospheric neutrino 
background for neutrinos energies E > 10 5 GeV. Typical results are shown in Fig. 9. These 
illustrate the possibility of observing AGN at various nadir angles by "neutrino astronomy" . 
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Table 1: The charged-current and neutral-current cross sections (in cm 2 ) for vN and vN 
interactions, where N = \{jp + n). 
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Figure Captions 



Fig. 1 Diagrammatic representation of the kx- factorization formula (0). At lowest order in 
as the gauge boson-gluon fusion processes, Vg — > qq, are given by the quark box shown 
(together with the crossed box). The variables k, k and k! denote the transverse momenta 
of the indicated virtual particles. 

Fig. 2 The total vN charged current cross section and its decomposition into components of 
different origin as a function of the laboratory neutrino energy. 

Fig. 3 The total vN neutral current cross section and its decomposition into components of 
different origin as a function of the laboratory neutrino energy. 

Fig. 4 The total vN cross section together with its charged current and neutral current compo- 
nents as a function of the laboratory neutrino energy. 

Fig. 5 The comparison of the total uN cross section calculated with and without the consistency 
constraint (|j) imposed. 

Fig. 6 The prediction for the total uN cross section obtained from a unified BFKL/DGLAP 
equation with (and, for comparison, without) the consistency condition imposed, com- 
pared to other recent calculations: |4| based on CTEQ parton distributions and || 
based on GRV dynamical partons |TB| . 



Fig. 7 A contour plot showing the x, Q 2 domain of the dominant contribution to the 
da / d ln(l/x)d log Q 2 for the total vN interaction at two values of the neutrino laboratory 
momentum: (a) E v = 10 6 GeV and (b) E v = 10 11 GeV. The 20 contours are such that 
they enclose a contribution of 5, 10, 15 • • • % of the above differential cross section. 



Fig. 8 The shadowing factor S of Q2"8] ) for two different initial neutrino fluxes incident at three 
different nadir angles on a detector. The angle 8 = 0° corresponds to penetration right 
through the Earth's diameter. The two curves on each plot show the shadowing factor 
with and without NC regeneration included. 

Fig. 9 The initial flux Iq(E) and the flux at the detector 1(E) for three different nadir angles cor- 
responding to three models for AGN neutrinos p3[ |24|, [25j . The background atmospheric 
neutrino flux is also shown. All the fluxes are given for muon neutrinos. The correspond- 



ing fluxes from [23, 24, 25 were given originally for muon neutrinos and anti-neutrinos, 



and their value has been divided by factor 2. 
Fig. 10 As for Fig. 9, but showing neutrino fluxes from gamma ray bursts PB1 and from a top- 



down model P7| .All the fluxes are given for muon neutrinos. The corresponding fluxes 



from [26, ^7| were given originally for muon neutrinos and anti-neutrinos, and their value 



has been divided by factor 2. 
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